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Abstract 

From extensive numeric diagonalizations of the SU{3) Perk-Schultz 
Hamiltonian with a special value of the anisotropy and different boun- 
dary conditions, we have observed simple regularities for a significant 
part of its eigenspectrum. In particular the ground state energy and 
nearby excitations belong to this part of the spectrum. Our sim- 
ple formulae describing these regularities remind, apart from some 
selection rules, the eigenspectrum of the free fermion model. Based 
on the numerical observations we formulate several conjectures. Us- 
ing explicit solutions of the associated nested Bethe-ansatz equations, 
guessed from an analysis of the functional equations of the model, we 
provide evidence for a part of them. 

1. Introduction 

Since the pioneering work of H. Bethe in 1931 the Bethe ansatz and its 
generalizations have proved to be quite efficient tools in the description of the 
eigenvectors of a huge variety of one-dimensional quantum chains and two- 
dimensional transfer matrices (see, e. g., |]1[ for reviews). Models with wave 
function given by this ansatz are considered exactly integrable. According 



to the Bethe ansatz the amphtudes of the wave functions are expressed in 
terms of a sum of plane waves whose wave numbers are given in terms of non- 
hnear and highly non trivial coupled equations known as the Bethe-ansatz 
equations (BAE). These equations in several cases, thanks some appropriate 
guessing on the topology of roots, are solvable in the thermodynamic limit 
providing the understanding of the large-distance physics. 

However the exact integrability is a property independent of the lattice 
size and the exact solution of the associated BAE for finite chains is an impor- 
tant step toward the complete mathematical and physical understanding of 
these models. Due to the complexity of the BAE up to our knowledge, only 
in two special cases some of the solutions are known analytically, namely, 
the trivial free-fermion case and the XXZ chain at the special anisotropy 
A = — 1/2 0, 1^. The solution in this last case is obtained by exploring the 
functional relations of the model |0. Even in the last case, although several 
exact properties of the wave function were conjectured 0] a complete and 
closed calculation of their amplitudes is still missing. In this paper we are 
going to present a new set of analytical solutions of BAE for finite chains. 
These solutions correspond to BAE of the anisotropic SU{3) Perk-Schultz 
model 0, or the anisotropic SU{3) Sutherland model 0, at a special value 
of the anisotropy. Contrary to the XXZ case the Bethe ansatz for this model 
is of nested type and the solution are going to be derived by generalizing the 
functional method originally applied to the XXZ chain. 

The paper is organized as follows. In section 2 we give the main defini- 
tions and formulate the corresponding BAE. In section 3 we state a set of 
conjectures that were obtained from extensive "experimental" work on exact 
bruteforce diagonalizations of the quantum chains. In section 4 we derive, for 
the Hamiltonian with periodic boundary conditions the functional relations 
and at a special value of the anisotropy some solutions for the eigenspectra 
are derived. In section 5 we present and test directly a set of solutions of the 
BAE, and explain partially the conjectures announced in section 3. Finally 
in section 6 we present our conclusions and a summary of our results. 

2. The SU{3) Perk-Schultz model 

The SU{3) Perk-Schultz model is the anisotropic version of the SU (3) 
Sutherland model with Hamiltonian, in a L-sites chain, given by 
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L-1 

HI = E + pHl,i (p = 0, 1), where (1) 

1 2 

Hi,j = - E E (EfE';'' + E^'^Ef - qE^Ef - l/g^f ) 

a=0 f)=a+l 

The 3x3 matrices E""^ have elements {E°'^)(.d = ^c^^ ^^"^ ? = exp(i?7) is the 
anisotropy of the modeL The cases of free and periodic boundary conditions 
are obtained by setting p = andp = 1 in (|lD, respectively. This Hamiltonian 
describe the dynamics of a system containing three classes of particles (0,1,2) 
with on-site hard-core exclusion. At g = 1 the model is S'?7(3) symmetric and 
for q ^ 1 the model has a [/(I) (g) ?7(1) symmetry due to the conservation 
of the number of particles of each specie. Consequently we can separate 
the Hilbert space into block disjoint sectors labeled by (no,ni,n2), where 
Hi = 0, 1, L is the number of particle of specie i (i=0,l,2). The Hamiltonian 
has a 5*3 symmetry due to its invariance under the permutation of distinct 
species, that implies that all the energies can be obtained from the sectors 
{nQ,ni,n2), where ^ '^i ^ "^2 and tiq + rii + n2 = L. Moreover in the 
special case of free boundaries {p = 0), the quantum chain (|I|) is also invariant 
under the additional quantum SU (3)q symmetry implying that all energies in 
the sector {n'(^,n[,n2) with Uq < n[ < n'2 are degenerated with the energies 
belonging to the sectors (77,0,^1,^2) with no < ni < 712, if ti'q < no and 
n'o + n[ <no + rii. 

The eigenenergies of the Hamiltonian (|l|) for p = or p = 1 in the sector 
{nQ,ni,n2) are given by 

m+ni 



J, ^ 1 I sm(%-W2) ^ sin(^, + r7/2) Y 

^ q sin(uj + 77/2) sin(Mj — 77/2)/ 



(2) 



where the variables {uj,j = 1, 2, no + ni} and the auxiliary variables 
{vj,j = 1,2, ...,no} are the roots of the coupled Bethe ansatz. These BAE 
are of nested type and in the case of periodic boundary they are given by 
(see e. g. 0, i) 

sin(Mfc + 77/2) 1^ = _"tt' sin(Mfc - Uj + r]) ^ sinjuk - vj - v/2) 
sm{uk - ri/2) \ fj^ sm{uk - Ui - r]) ^J-^ sm{uk - vj + ?7/2) ' 

smjvi -Vj + r]) "«+"^ smjvi - uj - rj/2) ^ _^ 
fj^ sm{vi - Vi - 7]) sm{vi - Uj + r]/2) 
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where k = 1, uq + rii and / = 1, no- 

In the case of free boundary the BAE are given by 

sin(Mfc + 77/2) 1 ^ sin(Mfc + Vi + r]/2) sm{uk -Vi + r]/2) 
sin(Mfc - 77/2) J f}^ sin(Mfe + f » - r]/2) sin{uk ~ Vi - t]/2) 
^o+rn sm{uk + Uj + T]) sin(-Ufc — Uj + rf) 
jJi j^k + Uj - ifj sin(Mfc - Uj - rf) ' 

no 



-r-r sin(i;; + Vi + rii) sm.{vi - Vj + r]) 
sin(t;; + Vi - r]) sm{vi - Vi - t]) 



no+ni 



-r-r sin(t;; + Uj + r]/2) sin(f; - Uj + 77/2) 
^-^^ sin(f; + - 77/2) sin(fz - Uj - 7//2) ' 

where = 1, tiq + 77,1 and / = 1, tiq. In the case of periodic boundaries 
the momentum P = ^ (/ = 0, 1, ...,L — 1) of the associated eigenstate is 
given by 

The solutions of the BAE are going to provide the eigenenergies of (|l|) if 
they correspond to non-zero norm Bethe states. The combinatory nature of 
the Bethe wave functions imply that solutions of (^) or (^) with coinciding 
roots produce null states. Nevertheless the requirement of non-coinciding 
roots does not necessarily ensure a genuine eigenvector, since conspicuous 
cancellation even in this case, can render a vector with null norm. Although 
all eigenenergies of the Hamiltonian can be obtained, apart from predicted 
degeneracies, by restricting to the sectors (tiq, rii, 71-2) with tiq < rii < n2, the 
Bethe ansatz implementation in its coordinate version is valid for arbitrary 
values of tiq, 7i2. However as we shall see in section 5, several solutions 
with non-coinciding roots for sectors out of the range Uq < rii < n2, but 
corresponding to null state can be obtained. In fact even in the XXZ chain, 
where the BAE are simpler, solutions with non-coinciding roots Q that corre- 
spond to null-norm states can be obtained when the number of roots n is out 
of the range n < L/2. In the case of the BAE for the XXZ chain, in a recent 



paper [|ri[] Baxter gives strong evidence that with suitable parametrization 



^We have found, for example, a continuous set of non-coinciding roots of the BAE 
for the periodic XXZ with L = 4 sites and n = 3. All this set of solutions give us zero 
eigenvectors (see also for further considerations). 
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the entire eigenspectra can be obtained from the non- coinciding roots of the 
associated BAE in the sector with the number of roots n < L/2. In a similar 
way we are going to assume on this paper that all distinct eigenenergies of 
([^) can be obtained from the solutions {fi}, with non-coinciding values 
in the sectors when uq <ni < n2. Conversely solutions in sectors out of this 
range, should necessarily be degenerated with energies occurring in sectors 
within the range, if they do not correspond to zero-norm states. 

3. Conjectures merged from numerical studies 

In this section we state a series of conjectures that are consistent with the 
exact brute-force diagonalization of the Hamiltonian (|l|) with free {p = 0) 
and periodic (p = 1) boundary condition at g = exp(2i7r/3). Some of these 
conjectures are going to be proved in the following sections. Let us consider 
separately the case of periodic and free boundaries. 

3a - Periodic chain. 

CONJECTURE 1. The Hamiltonian (|l]) with L sites at g = exp(2z7r/3) has 
eigenvectors (not all of them) with energy and momentum given by 

E, = -^(l + 2cos^), (6) 

with / being a subset of X unequal elements of the set {1,2, The 
number I has to be odd X = 2/c + 1 and the sector of appearance of the 
above levels is Sk = {k, k + 1, L - 2k - 1), < k < {L - l)/2. 

The lowest eigenenergy among the above conjectured values (^) is ob- 
tained for the particular set Iq'^^ = {1, 2, A;} U {L — k, ...,L}, since in this 
case all contributions — (1 + 2cos^) to (^) has the lowest possible values. 
The corresponding eigenstate has zero momentum and energy given by 

= - V (1 + 2 cos g^) = -2fc - 1 - 2^^"("^'f + ^^/"-^ . (8) 
."Tfc) L sm(7r/L) 

CONJECTURE 2. For arbitrary L = 3n + I {I = 1,2,3), the eigenenergy 
Eq"''' is the lowest one in the sector Sn, moreover if / 7^ 3, (L 7^ 3n) it is the 
ground-state energy of the model. 
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3b - Free boundaries. 

In order to state our conjectures let us define again the special set of 
sectors of the Hamiltonian with p = 0: 

S, = (Int(^), Int(^), L-k), k = 0,l, ...L. (9) 

Due to the quantum symmetry SU{3)q, distinct sectors show the same eigenen- 
ergies. For example, for L = 7 the sectors are 

5o = (0,0,7) 5i = (0,l,6) 

52 = (1,1,5) ^3 = (1,2,4) 

^4 = (2, 2, 3) 55 = (2,3,2) 

56 = (3,3,1) 57 = (3,4,0) 

and we have a special ordering 

5o C 5i C ^2 C ^3 C ^4 = ^5 D ^6 D ^7- (10) 

This means that, for example, all eigenvalues in sector 5*2 can also be found 
in sectors 5*3, 5*4 and 5*5, and on the other side all eigenvalues in sector S7 
also appear in sectors 5*6, S'5 and S4. Sectors 5*4 and 5*5 are totally equivalent. 
Let us call in this example the sectors Sq, Si, S2, S3, S4 as the LEFT sectors 
and S'5, 5*6, 5*7 as the RIGHT ones. This can be directly generalized to any 
L = 3n + 1 or L = 3r;, + 2, obtaining L — n left sectors and n + 1 right 
ones. In the case where L = 3n the sectors Sk with k = 0, l,...,2r;, — 1 
and k = 2n + 1, ...,L as the left and right sectors, respectively. The sector 
S2n = (^^5 n) is degenerated with two sectors S2n-i = (n — 1, n, n + 1) and 
S2n+i = (n, n + 1, n — 1) (S'2„_i = 5'2„+i) and can be considered either as a 
left or right sector. We state now the conjecture. 

CONJECTURE 3. Let L = 3n + / (Z = 0, 1,2). Then the Hamiltonian (HD 
with free boundaries (p = at g = exp(2z7r/3) has eigenvectors with energies 
given by 

^/ = -E(l + 2cos^), (11) 

with / is an arbitrary subset formed by k distinct elements of the set {1, 2, L— 
1}. Moreover if 5*^ is a left sector, then these eigenvalues appear in the sectors 
Sk, Sk+i, SL^n (^L-n+ifor/ = 0), and if 5*^ is a right sector, the eigenvalues 
appear in the sectors SL^n~i, Sl^u, Sk+i- 
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As a consequence of conjecture 3 the Hamiltonian (|T|) has the special 
eigenvalues 

and we are now in condition to formulate a remarkable conjecture. 
CONJECTURE 4. The lowest energy in the sector Sk is E^''^ or E^''"^^ 
depending if Sk is a left or a right sector respectively. 

The minimal value of E^''^ is obtained for k = L—n—1 and our "numerical 
experiments" induce the conjecture: 

CONJECTURE 5. The ground-state energy of the Hamiltonian (|^) with free 
boundary at g = exp(2z7r/3) is given by 

g„^i;a--)=2-L + n-^'°''''^''t^'/''-'. (13) 

sin(7r/2L) 

4. Functional relations for the anisotropic SU{3) Perk-Schultz 

model 

We are going to obtain analytically some of the conjectured results pre- 
sented in the previous section. Let us consider initially the periodic case when 
p = 1 in the Hamiltonian (0). The eigenenergies in the sectors with "par- 
ticle numbers" (no, rii, 112) are given by (H) where the Bethe roots {uj,j = 
1, 2, no + ni = 1712} and {vj,j = 1, 2, no = mi} are obtained by solving 
the BAE @. Bellow, to simplify the notation, we write A^^'' and A^^"* instead 
of Vj and uj, respectively. 

Defining the pair of sine-polynomials 

nil 

Qi{X) = llsm{X-Xf), 1 = 1,2, (14) 

the BAE (0) can be rewritten as 

Qi(Af + V) Q2{\f - r?/2) + Qi{\f - r,) Q^^xf + r,/2) = 

(j = l,2,...mi), (15) 
sin^(Al') + r//2) gi(Af + r//2) Q,{\t^ - ^) + 
+ sin^(Al') - r^/2) Qi{\f - Q^iXf + r^) = 

(A; = l,2,...m2). (16) 
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Since from the definitions Qi{Xf'') = for any Bethe roots = 1, 2), 
we should have the functional relations 

Qi(A + v) Q2{X - v/2) + Qi{X- v) Q2{X + = T2(A)Qi(A), (17) 
sin^(A + 7]/2) gi(A + 7]/2) Q2{\ - r^) + 

+ sin^(A - r//2) Qi(A - ^2) Q^iX + = Ti(A)Q2(A), (18) 

where T2(A) and Ti(A) are unknown sine-polynomials of the order m2 and 
L + mi, respectively. Shifting A— >A=F?7/2in (|T^ and inserting the result 
in 



18) we obtain 



sin^(A T r?/2) Q2(A ±7]) + sin^(A ± r//2) T2(A T V2)}Qi(A T r]/2) = 
sin^(A ± Ti/2) Q,{\ T 3^2) + Ti(A)}g2(A). (19) 

We now suppose that Qi{\ ± 77/2) and Q2(A) have no common roots, in this 
case: 

sin^(A T r//2)Q2(A ±77) + sin^(A ± r7/2)T2(A T v/^)} = 

= T±(A)g2(A), (20) 

sin^(A ± r]/2) Q^i\ t 3^2) + ri(A) = r±(A)gi(A T W2), (21) 

where T^(A) are sine-polynomialsQ of the degree L. The subtraction of equa- 
tions (pTD among themselves give us 

sin^(A + r//2) Qi{\ - 3r//2) - T+(A)Qi(A - r//2) 

+T-(A)Qi(A + r//2) - sin^(A - r//2) gi(A + 3r//2) = 0. (22) 



Similarly both equations (PD|) give us the relation 



sin^(A) sin^(A + r]) Q2{\ - 3r//2) 

- sin^(A + 7]) T-{\ - ri/2) Q^iX - ri/2) 
+ sin^(A - 7]) T+(A + r]/2) Q^^X + n/2) 

- sin^(A) sin^(A - r]) QiiX + 3ri/2) = 0. (23) 

Up to now our relations are valid for arbitrary values of the anisotropy 
1] and we now are going to restrict to the particular case 77 = 27r/3 (g = 

^ These polynomials are the eigenvalues of the transfer matrices corresponding to the 
fundamental representations of SU{3) in the auxiliary space. 
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exp(2i7r/3), where the several conjectures announced in Section 3 were ex- 
pected to be vahd. At this special value of the anisotropy we have the 
symmetry 

Qi{X-?,7]/2)=Qi{\-iT)=Qi{\ + 'R)=Qi{\ + ?,r]/2) / = 1,2, (24) 
and equations (^) and (pHf ) are given by 

0(A) gi(A - tt) - T+(A) gi(A - 7r/3) + T-(A) Qi{\ + 7r/3) = 0, (25) 

and 



-sin^(A)</)(A-7r)Q2(A 



TT 



sin^(A + 27r/3)T-(A - 7r/3) Q2(A - 7r/3) 



+ sin^(A - 27r/3)T+(A + 7r/3) g2(A + 7r/3) 
where 



0, 



(26) 



0(A) = sin^(A + 7r/3) - sin^(A - 7r/3). 



(27) 

The shifting A A ± 27r/3 in (p5| ) and (p6| ) show that these equations 
are equivalent to the linear matricial equations 



m 

(A^ 

-T+(A 



T-(A + f ^ 



-T+(A) 
0(A + f^ 



1^) T-(A 



27r 
3 



T~(A) 
-T+(A + f ) 
) 0(A-f) 



and 



0(A-7r 
-T+(A- 
T-(A-7r 



T-(A-f) 
-0(A-f) 

-T+(A - f ; 



-T+(A + f) 
T-(A + f) 
0(A + f) 



gi(A- 


vr) 


Qi(A- 


f) 


Qi(A + 


f) 


Q2(A- 


vr) 


Q2(A- 


f) 


Q2(A + 


f) 



0, (28) 



0. (29) 



respectively. In (^9[) we defined the new function Q2(A) = sin^(A) Q2(A). It 
is clear that T2(A + tt) = T2(A — tt) and consequently from (|20|) T±{X + tt) = 
T±{\ — tt). Equations ( pHf ) and (^) imply that non trivial solutions are 
obtained if the determinants of the matrices appearing in those equations 
vanish. Actually, by shifting A — A + vr in the determinant coming from (p9[) 



we clearly see that this last determinant vanishes if the one coming from (pSD 
also vanishes. 

The calculation of the general solutions T^(A) that render a null deter- 
minant is a quite difficult task, however simple solutions can be obtained 
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(rank 1) by imposing a proportionality between the columns of the matrix 
generating the determinant^, i.e., 

0(A) _ T-(A + 27r/3) _ -r+(A - 27r/3) 

(30) 



-T+(A) 0(A + 27r/3) T-(A-27r/3) ' 
-T+(A) _ 0(A + 27r/3) _ -T-(A-27r/3) 
T-(A) ~ -T+(A + 27r/3) ~ 0(A - 27r/3) 

We can verify the above relations are equivalent to the independent equa- 
tions 

T+(A) T-(A + 27r/3) = -0(A) 0(A + 27r/3), (31) 

T+(A) T+(A - 27r/3) = 0(A) T-(A - 27r/3). (32) 

In order to find solutions of these last equations, it will be useful to use the 
general relation 

L 

- 6^ = - ^■'^). ^ = exp(27rVi^) (33) 
i=i 

to write 

L 

0(A) = sin^(A + 7r/3) - sin^(A - 7r/3) = [] /z(A), (34) 

where 

/,(A) = sin(A + 7r/3)-c<;'sin(A-7r/3), (/ = 1,..,L). (35) 

Now consider any subset I of non-repeated integers of the set Jq = {1, 2, L}, 
and the complementary subset / , such that J J = Jq. We may try to solve 
(|31|)-(|33) by the ansatz 

^^(A) =to'^n/KA±27r/3) n /m(A), (36) 
ie/ me/ 

where are unknown constants. This ansatz imply 

T+(A) T-(A + 27r/3) = t+to H /^A + 27r/3) n lm{\) = 

1=1 m=l 

= t+to 0(A + 27r/3)0(A), (37) 



■^The idea to consider decreased rank in the functional relations was used previously in 
L2| to explain simple energy levels of a special case of the XXZ chain. 
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where ( 
straint 



was used in the last equality. The equation (^) imply the con- 

-1. (38) 



''0 ''O 



Also from (|3^) and (^ 



T+(A) T+(A - 27r/3) = {ttfm \{fi{\ + 27r/3) \{ /„(A - 27r/3) = 
= (t+)V(A)T-(A-27r/3)/to, (39) 



and (|32D imply, by using (^8]), that 



(^0 



1, U 



(40) 



Then the ansatz ( |35D with (|38| ) gives us a set of solutions for T^(A), that when 
inserted in the matricial equations (^8]) and (^) will provide the function 
Qi(A) and Q2{^)- The zeros of these last functions are the Bethe-ansatz 
roots and the eigenenergies are calculated by using in (Q) the roots of Q2(A). 
Instead of calculating the energies through this procedure, we are going to 
calculate them using the transfer matrix eigenvalues T~(A). From (|^) and 
the definition of Q2{^) it is not difficult to obtain the relation 



E 



V3 d / Q2(A) 



In 



2 dX \Q2{-X) 



A=7r/3 



(41) 



On the other hand let us expand (p3D with ?7 = 27r/3forA = ?7 + e e-Cl. 
The terms of the lowest order give us the relation 



d 



Q2(A) 

dX "Ag2(-A) 



In 



A=7r/3 



L 

71 



d 

dX 



lnr-(A)| 



A=7r/3; 



(42) 



that from (|41|) provide the simple result 



(43) 



The eigenenergies associated to our solutions T (A) are then obtained by in- 
serting (p6D in (1^, and we obtain after some simple algebraic manipulations 



E 



-L+^_(l + 2cos(-^)) 

m£7 



2cos( — )), 



(44) 
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where we used the formula 

L 

^ (1 + 2 cos(27r//L)) = ^(1 + 2 cos(27r//L)) = L. (45) 

Also the zero-order term in the same expansion of (p3D give us 
T-(7r/3) g2(vr/3) sin^ - r//2) 



sm 



(27r/3) g2(-7r/3) sin^ + ^2) 



n , = exp(.P), (46) 



where from (1) P is the momentum associated to our solution T (A) given 
in (pBf). Inserting (pB) into (EUf) we obtain after some simple calculations 



exp(^P) = (-l)^+it-exp(-^ 5:_m) = exp(^ J]/), (to)^ = 1. 

(47) 

5. Analytic solutions of the Bethe ansatz equations 

As we discussed in the previous section, at least in the periodic case, the 
Bethe ansatz roots, corresponding to the eigenenergies (44) we observed, can 
be obtained from the expansion of Q2(A) given in (0), derived by solving 
([29| ) with T^(A) given by the ansatz (^). Distinctly in this section we are 
going to present in a direct way a set of guessed solutions {ui.Vj} of the 
BAE that gives the energies conjectured in section 3. We show that they are 
correct by a direct substitution into the BAE. We present solutions of the 
BAE for the periodic and free boundaries cases. As we conjectured in section 
3, in the case of periodic boundaries there exist some selection rules in the 
spectrum composition (see conjecture 1). At the end of this section we are 
going to explain partially this conjecture. 

Let us consider separately the periodic and free boundary case. 
5a. Periodic case. 

The BAE (|^) at ?7 = 27r/3, expressed in terms of the variables 

sin(Mfc-7r/3) sin(t;; - 7r/3) , 

= ~f \ — 7^' = ~f , 7^' 

sm^Mfc + 7r/3) sm(t>; + 7r/3) 

with k = 1,2, ...,no + rii and / = 1, 2, rio are given by 

/ii 1 + yr + Jii i + Xk' + VjXk' 

(j = l,2,...no) 
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and 

/ii 1 + Vj' + ^kVj' Jii i + Xk' + XkXk' ^ 

(/c = 1,2, . . .no + ni). 

Let us fix 2no + rii = L. Our guessed solutions are obtained by considering 
{x/i, [k = 1, . . . ,no + ni,l = 1, . . . , uq) as an arbitrary permutation of 
{Uju"^, ...,uj^}, where uj = exp{2ni/L). In this case, the left side of equation 
(|49| ) takes the form: 



L 



.^^L+l-Q - ' ' (51) 



1 + + 



as 

\L I / i\L+l„,L 



,L+i(i + %)" + (-ir+X 



Using the identity (|33D and the fact that yj=lwe can rewrite this product 

V l + (_l)L+i(i + ^^.)L 1 • (52) 

It is evident that the second BAE is also satisfied due do equality = 1. 

Consequently we have found a subclass of solutions for the nested Bethe 
ansatz equations. These solutions are characterized by the subset / with 
unequal elements of the set Iq = {1, 2, L}, and have the energy 

no+rii 

Ej = - (1 + Xfc + 1) = - ^(1 + 2 cos(27rZ/L)) (53) 

k=l l£l 

and momentum 

no+ni 1 r) 

Pj= -Hx,) = -Yl- (54) 
k=i ^ ^ lei 

Comparing the above relations with relations (P) and (^ we verify that our 
guessed solutions are consistent with the conjecture 1. It is not clear however 
if the corresponding wave function is not a zero vector. 
5b. Free boundary case. 

The BAE (^ at r/ = 27r/3, expressed in terms of the same variables and 
yi with k = 1,2, ...,no + rii and / = 1, 2, no are given by 

1 + 1/i + yjVj' \fyj + vr + yjVj 



no 

V 1 + yy + y^yj' J\ 1 + + 



13 



\l + Xk' + ViXk' J\ l+Vi + Xk' J 



k'=l 



(j = l,2,...no) 



and 



^ / 1 + Xfc + XkVj' \ f Xk + Vf + XkVj' 
/ii ^ 1 + Vj' + ^kVj' J\ 1 + Xk + Vf 

V 1 + Xfc/ + XfcXfc' /V 1 + Xk + Xk' J ^ 
{k = 1,2, ...rio + ni). 

Now let us fix 2nQ + rii = L — 1. Our guessed solutions are now given 
by the set {xk, yi}{k = 1, . . . ,no + riij = 1, . . . , no) formed by an arbitrary 
permutation of {uj,uj'^, ...,uj^^^}, where u = exp(i7r/L). Using the identity 



ia + uj"^)il/a + ujn ^ b^-' jb' - 1) ja'^ - 1) 
iii (6 + cu™)(l/fe + c^-) a^-i (a2 - 1) _ 1) ' ^ ^ 



and the fact that = 1 we can easily verify that the BAE ( |55| ) and ( |56D are 
satisfied. 

As in the periodic case, we have found a subclass of solutions for the 
nested BAE. These solutions are characterized by a subset /c{l,2,...,L— 1} 
and have the energy 

no+ni 

Ei = - J2 (l + Xfc + x^^) = -5](l + 2cos(7r//L)). (58) 

fc=l /£/ 

All these solutions which are consistent with conjecture 3, so we think that 
corresponding Bethe wave function is not a zero vector. 

Finally in order to conclude this section we intend to explain partially the 
selection rules formulated in Conjecture 1 for the periodic case. We are going 
to this by exploiting our solutions (|36|) for T'^iX) of the functional relations 



of section 4 with the help of some ideas developed in the papers HlO 



Inserting our solutions ( p6D for T'^iX) into equation (|25| ) we obtain 

n fiw Qi(^ -^)-Ki[fi(^+ V3) n fmW - 71/3) + 

1=1 lei rnel 

to n //(A - 27r/3) n fmiX) + vr/3) = 0. (59) 

lel m&I 
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Dividing by the common factor Yim^i fmW we obtain 

i^i(A)gi(A-7r)+f]Fi(A+27r/3)gi(A-7r/3)+l]^Fi(A-27r/3)gi(A+7r/3) = 0, 

(60) 

where = — (^^ = 1) and 

F^iX) = l[MX). (61) 

lei 

On the other hand the solution ( |36D for T^{X) brings ( p5D into a similar 
functional equation: 

i^2(A)Q2(A-7r)+f]2F2(A + 27r/3)g2(A-7r/3)+l]F2(A-27r/3)g2(A+7r/3) = 0, 

(62) 

where 

F2(A) = sin^An/m(A). (63) 

Let us consider the case where L ^ 3n. In this case since P = 
{j = 0, . . . , L — 1), equation ( ^Tj) gives = 1, and consequently f2 = 1 in 
M) and (1621). 



We intend to argue now that there exist pairs {gi(A), g2(A)}, satisfying 



([60| ) and ( |62D with = 1 which lead to "physical" solutions for the nested 
BAE d^), i. e., solutions which are inside the usual bounds uq < rii < n2 or 
equivalently 

deggi <degg2-deggi <L-degg2. (64) 

First of all, we have special solutions for (^) and ( p^ which can be written 

as 



gi(A) = QispecW = n /-(A + vr), g2(A) = Q2spec{X) =UM^ + '^) (65) 

m&I lei 



Let us check these formulas inserting them into equations (|60|) and (|62D. 
The left side of equation (|^) becomes (see (|3^) ) 



n/KA) n/-(^) +n/KA + 27r/3) n/-(A + 27r/3) + 
lei mei lei mei 
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+ n/KA-27r/3) n/-(^-w3) = 

l&I m&I 

= sm^(A + 7r/3) - sin^(A - 7r/3) + sin^(A + tt) - 
- sin^(A + 7r/3) + sin^(A - 7r/3) - sin^(A - tt) = 0. 



Similarly the left side of equation (p^) becomes 



sin^(A) (sin^(A + n/3) - sin^(A - ir/S)) + 

+ sin^(A + 27r/3) (sin^(A + vr) - sin^(A + 7r/3)) + 

+ sin^(A - 2n/3) (sin^(A - n/3) - sin^(A - tt)) = 0. 

Let < X < L is the number of elements of the set I. Then degrees of 
these special solutions Qi and Q2 are L —X and X respectively. Inequalities 
(|6^ ) for these pairs become 

L-I<2I-L<L-X, (66) 

which is the same as the equality 2L = 31. It is not enough for our purposes, 
especially for L 7^ 3n, so we have to look for additional solutions. They do 
exist due to the fact that the matrices in equations p8| ) and ( ]29|) for Qi and 
Q2 has rank 1. 

According to the analysis of functional equations of type (pOD or ( ]62[ ) 



made in some previous papers ||T0[ it was noticed that equations of this type 
have some conjectured properties that we are going to accept. If in ( pDj ) or 
(|6^ ) Fi{\) [i = 1,2) have a trigonometric form Fi{\) = n^isin(A — a^). 



of degree N, in general there exists a trigonometric solution of the form 
QiW = rijli sin(A — 6j) of degree m. This degree depends on the value of 
Q appearing in the equation. In particular if f2 = 1 then m = N/2 + 1 for 
even and m = {N — l)/2 for odd. Only for special choices of Fi{X) 
this degree can be decreased. We call these solutions Qigen, Q2gen as general 
ones. 

Due to dll) and (H) we have deg Fi(A) = X and deg F2(A) = 2L-I. 



If we chose X even then 2L — X is also even and the equations ( IBOD and 
(|6^ ) have trigonometric solutions for Qi and Q2, with degQi = X/2 + 1 and 
degQ2 = {2L — X)/2 + 1. On the other hand for odd values of X we have 
degQi = (X - l)/2 and degQs = {21-1- l)/2. 

Before consider arbitrary values of L let us restrict initially to the particu- 
lar case L = 7. In table 1 we list the predicted degrees of the sine-polynomials 
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Qi and Q2- We underline pairs Qi, Q2 which satisfy the inequahties ( |51D and 
in the last column of this table we present the eigensectors where we expect 
to find the predicted simple energy levels. 



X 


deg Qigen 


deg Q2spec 


deg Qispec 


deg Q2gen 


sector 





1 





7 


8 




1 





1 


6 


6 


(0,1,6) 


2 


2 


2 


5 


7 




3 


1 


3 


4 


5 


(1,2,4) 


4 


3 


4 


3 


6 




5 


2 


5 


2 


4 


(2,2,3) 


6 


4 


6 


1 


5 




7 


3 


7 





3 


(0,3,4) 



First of all we see that only odd X leads to "physical" solution. This fact 
is consistent with the results of our "experimental" observations formulated 
in Conjecture 1. 

We see further that for small X the "physical" solution is a pair consisting 
of a general solution Qigen and a special one Q2spec- For odd X = 2fc + 1 we 
have deg Qigen = (^ ~ l)/2 = k and deg Q2spec = T = 2k + 1. Inserting 
these formulas into inequalities (0) we get k<k + l<L — 2k — 1. For 
L = 3n + I (Z = 1, 2) we obtain the upper boundary for k: 

k<n + (67) 
3 

On the other side for X large enough we combine a special solution Qupec, 
which has degree L—X= L — 2k — l and a general one Q2gen which has degree 
(2L-X-l)/2 = L-k-l. Inequalities (^ become L - 2fc - 1 <k<k + l 
Taking L = 3n + / (/ = 1, 2) we obtain now the lower boundary for k: 

k>n + (68) 

3 

There is no holes between ( p7| ) and ( pBD so we have "physical" solution 
for every odd X and corresponding energy level have to be in sector (/c, k + 
1, L — 2fc — 1). This explains Conjecture 1! 

The case L = 3n is more complicated and we did not derive similar 
results. 
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6. Summary and Conclusions 



Although the exact integrability is a property independent of the lattice 
size, the exact solution of the associated BAE for finite chains were know 
in very few cases. The XXZ chain at the special value of the anisotropy 
A = (g + g~^)/2, q = cxp(i27r/3) is one of these examples. Motivated by this 
result we made extensive numerical calculations for the SU (3) generalization 
of the XXZ chain, namely the SU (3) Perk-Schultz model, also at the spe- 
cial anisotropy q = exp(i27r/3). Surprisingly, as we stated in section 3, the 
numerical results reveal that many of the eigenenergies (not all of them) are 
expressed as combinations of simple cosines and, apart from some selection 
rules, are quite similar to the energies of a free fermion chain (or XXZ at 
A = 0). 

Our numerical results indicate the five conjectures presented in section 
3. The first two conjectures concerns with the periodic quantum chain and 
gives the exact expression for the energy and momentum of several eigen- 
functions. In several sectors the lowest energy value is also predicted. In 
order to explain these results analytically we present in section 5 a set of 
BAE solutions that are consistent with the conjectured energies. However 
the set of solutions we obtained is larger then the conjectured one. This 
imply that some of our solutions, although having non-coinciding roots are 
unphysical, corresponding to a zero vector, since the associated energy is 
missing in the eigcnspcctrum. These missing BAE solutions appear in the 
sectors (no, ni, 77,2) not satisfying the bound Uq < rii < n2 < 2L/3. From the 
functional relations derived in section 4 we were able to explain at least for 
the cases L ^ 3n the selection rules appearing in the conjecture 1. In the case 
L = 3n the degree of the trigonometric solutions of the functional equations 
are more difficult to predict and we could not explain the conjecture 1. 

The last three conjectures concerns the eigenspectra of the Hamiltonian 
with the quantum symmetry SU{3)q, i. e., free boundary case. These con- 
jectures shows no selection rules in contrast with the periodic case. Again 
in this case we present a set of solutions of the BAE sharing the same ener- 
gies as those of conjecture 3. The functional relations in this case are more 
complicated and we leave this analysis for a future work. 

Finally it is interesting to mention that the finite-size corrections obtained 
from the conjectured eigenenergies of section 3 give us some conformal di- 
mensions of the underlying conformal field theory (CFT) governing the large- 
distance physics of the model. As a consequence of the conformal invariance 
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of the infinite system these eigenenergies should behave as 



E = eooL + ^Vs{12xo-c)+oiL-^), (69) 

in the periodic case, and 

E = e^L + f, + ^Vs{24x: - c) + o{L-^), (70) 

in the open boundary cases. In the above expression Coo and fs are the 
energy per site and surface energy in the bulk limit, Vs is the sound velocity, 
c is the central charge and the conformal dimensions governing the 

power-law decay of correlations in the periodic and open chain cases. 

In the periodic case, the conjecture 2 gives the asymptotic behaviour for 
the lowest eigenenergies 

E = e^L-^Vs{~2) + 0{L-'') for L = 3n, (71) 

E = e^L-^Vsl + 0{L~'') for L ^ 3n, (72) 
oL 3 

where eoo = ~(| + "v) ^"^^ ^■^ ~ '"^^ inferred from the lowest eigenen- 
ergy with momentum ^ of conjecture 1. 

The underlying f/(l) ® f^(l) CFT governing these quantum chains is 
expected to have a central charge c = 2 and when formulated in the torus 



should have the conformal dimensions 



a;(ni, 712; "^i, "I2) = Xp{n\ - nin2 + n\) + 7777-(m'i + mim2 + m\), (73) 



12a;p 



where Xp is related with the compactification ratio and ni,n2,mi,m2 are 
expected to be integers. Assuming c = 2 in ( p5| ) and comparing with relations 
[1|) and (^) we obtain for the predicted lowest eigenenergies, the associated 



dimensions x = | for L = 3n and x = | for L 7^ 3n. From (^) these 
dimensions can be identified with x(l, 1; 0, 0) = Xp = | and x(|, — |; 0, 0) = 
^ = by taking Xp = ^ = ^. The fractional values in the last case happens 
because the ground state for lattices with sizes L ^ 3n does not represent, 
in the bulk limit, the true vacuum of the CFT, since it contains topological 
defects. 
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In the case of open boundaries conjecture 5 give us for the ground state 



E = e^L + f,- ^v,{-2) + 0(L-3), (74) 

where Coo and Vg was aheady obtained in the periodic case and fs = 3/2. 
Comparing ([7^ ) with ( |70D we obtain c = —2. This can be understood since 
the quantum chain with open boundaries is SU{3)q symmetric, with q = e*'', 



Tj = and the expected |]T5| conformal anomaly in this case is c = 2- 
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m(m+l) ' 



where m = = 2. Similar analysis can also be done for the excited states. 
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TABLE CAPTION 

Table 1 - Degrees of the polynomials Qi and Q2 coming from the possible 
solutions for L + 7. The special solutions Qupec and Q2spec are given by (pSD 
and the general ones Qigen and Q2gen are discussed in the text. 
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